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Abstract 

Employing a limiting case of a conjecture for constructing piecewise 
separable-variables functions, the elements of the Pseudoanalytic Func- 
tion Theory are used for numerically approaching solutions of the for- 
ward Dirichlet boundary value problem, corresponding to the Electrical 
Impedance Equation in the plane, when the electrical conductivity is an 
arbitrary non- vanishing function, fully defined within a bounded domain. 
The new method is studied considering a variety of examples when the 
bounded domain coincides with the unit circle, and it is also included a de- 
scription of its behaviour in non-smooth domains, selecting special cases 
that do not require additional regularization techniques, for warranting 
the convergence of the approach at the non-smooth regions, when certain 
requirements are fulfilled. 

1 Introduction. 

The study of the Electrical Impedance Equation 

V • (crVw) = 0, (1) 

where a represents the conductivity function and u denotes the electric poten- 
tial, is fundamental for the proper understanding of a wide variety of physical 
problems, among which we find the Electrical Impedance Tomography, correctly 
stated in mathematical form by A. P. Calderon |5] in 1980. Indeed, most of the 
algorithms that approach solutions for this inverse problem, within a bounded 



domain $7 (R^) , are based on iterative methods that examine solutions for the 
direct problem, and introduce a certain kind of variations in the conductivity u, 
attempting to minimize the difference between the approached solution u, and 
the boundary condition Mc|r(see e.g., the classic work [14,). Here F denotes the 
boundary of the domain Q. 

Yet, specifically talking about the Electrical Impedance Tomography in the 
plane, the mathematical complexity of ([I]) has posed so strong challenges, that 
the problem is still considered ill posed. This is, when small variations of the 
conductivity a are introduced for minimizing the error, most of the classical 
numerical methods do not achieve to decrease such error under a certain limit. 

At this point, it is convenient to empathize that the majority of these 
methods, are based upon variations of the Finite Element Method, one of the 
finest tools for solving partial differential equations in bounded domains (M^) . 
Therefore, the usage of techniques upcoming from completely different branches 
of the Applied Mathematics, could well show up new information about the be- 
haviour of the solutions for the forward Dirichlet boundary value problem of 
([T]) in the plane, which eventually could allow to propose new techniques for 
analysing its inverse problem. If such becomes true, its common classification 
of ill posed problem could be reconsidered. 

In this sense, the discovering of the relation between the two-dimensional 
case of ([T]) and the Vekua equation [T3] , independently achieved by V. Kravchenko 
in 2005 [S], and K. Astala and L. Paivarinta in 2006 P], possesses special rel- 
evance, since it opened a new path for analysing the solutions of ([!]) in the 
plane, from the point of view of the so-called Taylor series in formal powers [2]. 
As a matter of fact, the Vekua equation had been deeply studied in a variety 
of interesting works, published almost five decades before its relation with ([!]) 
would have been first noticed. Indeed, two of the most important and complete 
works about the Vekua equation, were published by L. Bers [2] in 1953, and I. 
Vekua [l3] itself, in 1962. 

The number of publications that came after discovering the cited relation is 
long, and they are all interesting (see e.g. [B] and [7]). Yet, we shall remark 
that only some of them can be directly employed in Physics (see e.g. [S]), since 
it is not always clear how to adapt the elements of the Modern Pseudoanalytic 
Function Theory to the experimental physical requirements. 

The current pages intend to make a positive contribution in this direction. 
We will not reach the study of the Electrical Impedance Tomography problem, 
but we will provide the elements that allow the employment of the Modern Pseu- 
doanalytic Function Theory, in a variety of cases that can be easily identified 
with physical experimental models. 

More precisely, after providing the necessary elements of the pseudoanalytic 
functions, we expose the numerical method that allows the construction of the 
formal powers, which eventually, will provide a complete set of solutions for 
solving the forward Dirichlet boundary value problem of ([T]), when Q coincides 
with the unit circle. The performance of the numerical method is tested by em- 
ploying one example of conductivity that possesses a separable-variables form, 
and for which an exact solution is known. Indeed, that exact solution will be 
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imposed as the boundary condition Mc|r, so we can estimate the accuracy of the 
approaching by a measure £, that coincides with the Lebesgue integral over the 
boundary T . 

Subsequently, we use a Conjecture |T2] for constructing piecewise separable- 
variables conductivity functions, a requisite for studying the solutions of ([T]) from 
the point of view of the Pseudoanalytic Function Theory |5] . We examine the 
effectiveness of this idea solving once more the boundary value problem posed 
previously, but employing the piecewise conductivity in lieu of the original cr, 
making a comparison between the errors £ reached in both cases. 

Immediately after, we present a Proposition stating that any conductivity 
function a : $7 (M^) M, can be considered the limiting case of a piecewise 
separable-variables function, at every point (a;,y) € O. As a matter of fact, this 
Proposition is the vertebral column of this work, because it allows the study of 
the forward Dirichlet boundary value problem of ([T]), for the cases when the ex- 
act representation of a is known, but does not posses a separable- variables form. 
Moreover, for the first time in the literature dedicated to the Applied Pseudo- 
analytic Function Theory, the numerical method will be employed for analysing 
conductivities upcoming from geometrical distributions, a very important fact 
for physical applications. 

We close this work with a review of the behaviour of the method, when 
the domain does not coincide with the unit circle, and possesses avoidable 
discontinuities in the derivative of the parametric curve describing P. To better 
analyse this case, we selected three kind of conductivities, two of them rising 
from geometrical distributions, that do not require any additional regularization 
technique, in order to warrant the convergence of the approached solutions in 
the corner points, when some certain conditions are fulfilled. 

Notice that, even a formal comparison of this new technique with an ade- 
quate variation of the Finite Element Method is in order, on behalf of briefness 
we refer the reader to the results reported in [7], where an accurate likening 
of this type was performed. The reader will find that, nonetheless the contrast 
was made employing separable- variables conductivities, the results can be easily 
extended to the cases treated in the current pages. 

2 Preliminaries. 

Following pi, let the complex- valued functions F and G satisfy the condition 

Im (FG) > 0, (2) 

where F denotes the complex conjugate of F: F = ReF — ilvtvF, and i is the 
standard imaginary unit: = —1. Thus, any complex- valued function W can 
be expressed by means of the linear combination of F and G: 

W ^(l)F + tPG, 

where and "0 are purely real functions. Two complex functions that fulfil 
Q shall be called an (F, G)- generating pair. Bers |2j introduced the (F, G)- 
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derivative of the function W according to the expression: 

5(i.,G)W^= (5.0)F+(9,^)G. (3) 
This derivative will exist if and only if 

(chcb) F + G = 0, (4) 

where 

dz^dx-idy, ch = dx+idy. 

Notice that these operators are classically introduced with the factor ^ , but 
it will result more convenient to omit it in this work. 
Introducing the functions 



A 



Fd.G-Gd^F Fd^G- Gd^F 

(F,G) - 7^^' "(^'G) 



FG-GF ' ' ' FG-GF 

_Fd,G-Gd,F _ GchF-Fd^G ... 



the expression of the {F, G)-derivative ([s]) will turn into 

d(F.G)W^d,W-A(^F,G)W-B^F^G)W, (6) 

and the condition Q will be written as 

d^W-a(^F^G)W^b^F,G)W = 0. (7) 

The functions defined in ([5]) are called the characteristic coefficients of the 
generating pair {F,G), whereas the functions W, solutions of the equation ([t]), 
are named (F, G)-pseudoanalytic functions. As a matter of fact, the equation 
^ is known as the Vekua equation 13J, and it is the foundation of the present 
work in many senses. 

The following statements were originally presented in [2] and [8] . They have 
been adapted here for the purposes of this work. 

Theorem 1 The elements of the generating pair (F, G) are {F, G) -pseudoanalytic: 

d{F.,G)F = d{F,G)G = 0. 

Remark 1 Let p be a non-vanishing function within a bounded domain (M^) . 
The functions 

Fo=P, Go = -, (8) 
P 

constitute a generating pair, whose characteristic coefficients are 

^(Fo,Go) = 0(^0, Go) = 0; 

BiFo,Go)=P2^9,p, (9) 

biFo.Go) =P~^Ckp. 
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Definition 1 Let {Fq, Gq) and {Fi, Gi) be two generating pairs of the form 
and let their characteristic coefficients satisfy the relation 

B(Fi.,Gi) = -bi^F^^Go)- (10) 

Thus, the pair {Fi,Gi) will be called a successor of the pair (Fq,Gq), whereas 
{Fq, Gq) will be called a predecessor of [Fi, Gi). 

Definition 2 Let 

{{Frn,G„,)}, m = 0,±l,±2,... 

be a set of generating pairs, where every {Fm+i, Gm+i) is a successor of{Fm, Gm) 
Therefore, the set {{Fm,Gm)} will be called a generating sequence. Moreover, 
if there exist a number c such that {F„i,Gm) — {Fm+a Gm+c) the generating 
sequence will be periodic, with period c. 

Finally, if {F,G) — {Fq,Gq), we will say that the generating pair {F,G) is 
embedded into the generating sequence {{Fm,Gm)}- 

Theorem 2 Let (F, G) be a generating pair of the form ^§^, and let p be a 
separable-variables function: 

p = Pi{x)p2{y), 

where x,y G M. Thus {F,G) will be embedded into a periodic generating se- 
quence, with period c — 2, such that 

P _ P2{y) ^ _ . Pijx) _ 

Pi{x) P2[y) 

when the subindex m is an even number, and 

1 

Fm = Pi {x)p2 (y) , = — ^ ; 

Pi(x)P2{y) 

when m is odd. 

Furthermore, if particularly pi{x) — 1, it is easy to verify that the generating 
sequence in which {F, G) is embedded will be periodic, but with period c = 1. 

L. Bers also introduced the concept of the (Fq, G'o)-integral of a complex- 
valued function W. We refer the reader to the specialized literature [5] and [S] 
for a detailed description of the necessary conditions for its existence. In the 
current pages, every complex function contained into an (i^Oj Go)-integral will 
be, by definition, integrable. 

Definition 3 Let {Fq^Gq) be a generating pair of the form Q). Its adjoin 
generating pair {Fq,Gq) is defined according to the formulas 

F*=-tFo, G* = -iGo. 



5 



Definition 4 The (Fq^Gq) -integral of a complex-valued function W (when it 
exists l^) is defined as: 

I Wd(Fo,Go)Z = FfiRe I G*Wdz + GaRe f F*Wdz, 

where rj is a rectifiable curve within a domain in the complex plane. Specifically, 
if we consider the {F^f^Go) -integral of the {Fq, Gq)- derivative ofW, we will have 
that: 

9(fo,Go)W^d(Fo,Go)^ = -^{zo)F{z) - ^{zo)G{z) + W{z), (11) 

where z = x -\~ iy, and zq is a fixed point in the complex plane. According to the 
Theorem^ the {Fq, Go)-derivatives of Fq and Gq vanish identically, thus the 
expression fTll) can be considered the {FQ,GQ)-antiderivative of d^p^ Q^-^W . 



2.1 Formal Powers. 

Definition 5 The formal power {oq, zq] z) belonging to the generating pair 
{Fm,Gm), with formal degree 0, complex constant coefficient uq, center at zq, 
and depending upon z = x -\- iy, is defined according to the expression: 

Z^^ (ao, zq;z)^ XF„^{z) + fiG,n{z), (12) 

where A and fj, are complex constants fulfilling the condition: 

XFraizQ) + ^Gto(zo) = Uq. 

The formal powers with higher degrees are approached according to the re- 
cursive formulas: 

("-n, Zq] z) ^ U ^^Ji^^ (a„, Zq] z) d(F^^G^)Z, (13) 



where n ~ 1,2,3. .... Notice the integral operators in the right-hand side of the 
last expression are {Fm,Gm,)-o.ntiderivatives. 

Tiieorem 3 The formal powers posses the following properties: 

1. Every zln' (a„, zq; z) , n = Q, 1,2, ... is an {Fm, Gm) -pseudoanalytic func- 
tion. 

2. Let a„ = a'^ + m", where a'^,a'^ CE M. The following relation holds 

Zi") (a„, ZQ- z) = (1, zq- ) + ^Z^,"' {i, zq- z) . (14) 

3. Finally 

lim (a„, Zq; z) = a„(z - zq)"- (15) 
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Theorem 4 Let W he an {F„nGm)-pseudoanalytic function. Then it can be 
expressed in terms of the so-called Taylor series in formal powers: 

oo 

W = J2zi:Han,zo;z). (16) 

n=0 

Furthermore, since any (Fm, Gm)-pseudoanalytic function W accepts this expan- 
sion, (16) is an analytical representation of the general solution for the Vekua 
equation 

2.2 The two-dimensional Electrical Impedance Equation. 

Let us consider the equation ([ij in the plane: 

V • (crVw) = 0. 

As it has been shown in several previous works (see e.g. f8| and [TT|), if cr can 
be expressed by means of a separable-variables function: 

cr{x,y) = ai{x)a2{y), 

by introducing the notations 



W = — iy/adyU, 

P = (VoT) Vo^; 



(17) 



the equation ([!]) will turn into the Vekua equation 



for which the functions 



d^W - —W = 0, (18) 



Fo=P, Go = -, (19) 
P 



conform a generating pair. 

From (17 1, and according to the Theorem [2] it is possible to verify that 
this pair is embedded into a generating sequence, with period c = 2, for p is a 
separable-variables function. 

2.3 A complete set for the Dirichlet boundary value prob- 
lem of the two-dimensional Electrical Impedance Equa- 
tion. 

An explicit generating sequence allows the construction of the formal powers 
( 13 ), and in consequence, the approaching of the general solution of ([T]), accord- 
ing to the relations ([Tt]). 

Indeed, a very important relation between the solutions of ([!]) and of ( 18 1 
was elegantly posed in [6] , and this relation will play a central role in the present 
work. 
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Theorem 5 Let us consider the set of formal powers 



{z(")(l,0;z), Z(")(i,0;z)r , 

^ J n—Q 



corresponding to the generating pair (19), and let 51 (]R^) he a bounded domain, 
with boundary T, such that G 51 but ^ F. Then the set of functions defined 
on T: ^ 

|i?e4"^ (l,0;z)|r, i?eZ^"^ (i, 0; z) |r|" , (20) 

conforms a complete system for approaching solutions of the forward Dirichlet 
boundary value problem of |7p. 

In other words, according to the second property of the Theorem |3] when a 
separable-variables conductivity function a is given within a bounded domain 
ft, and a boundary condition Mc|r is imposed for the solution of (Tj), it will be 
always possible to construct a finite set of functions, subset of (20), such that 



/ (^"clr -X]<Re^^"^ (1,0; z) |r + Re^"^ {i,0;z) |r^ 



dl < £, 



(21) 



where £ > and / G F. 



2.4 Construction of a piecewise separable-variables func- 
tion. 

Approaching solutions of the forward Dirichlet boundary value problem for ([T]), 
by employing formal powers, has proven its effectiveness in a variety of works 
(see e.g. [S] and [?])• Yet, how to apply those methodologies when the conduc- 
tivity (J is not represented as a separable-variables function, remains an open 
question. 

A possibility for studying these cases could be to introduce an interpolating 
method that, given a set of conductivity values defined into a bounded domain in 
the plane, it can be able to approach a piecewise separable-variables function. 
One of the first proposals in this direction was posed in [TT], and the next 
paragraphs will show that, even it is a basic idea, it can well serve to our main 
objectives. 

Consider a bounded domain f2 (a unitary disk, for instance), and divide it 
into a finite number of subsections, taking care that the center zq of the formal 
powers (see the Definition [5| does not reside onto the boundary of two or more 
subsections. On behalf of simplicity, let us make the division by employing a 
finite set of parallel lines to the y-axis, equidistantly located one to each other, 
and let us fix zq = 0. 



A simplified illustration of this steps is plotted in Figure 1(a) For this 
example, we consider K — 3 subsections, delimited by the set of if + 1 y-axis 
parallel lines 

{X(o),X(i),X(2),X(3)} . 
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(a) A circular sectioned do- (b) Points on the lines 
main. on which the interpolat- 

ing functions will be 

constructed. 

Figure 1: Simplified illustrations of the steps for approaching a piecewise 
separable- variables conductivity function. 



Supposing that the values of the electrical conductivity are defined at every 
point inside the domain Q, let us trace a straight line within every subsection, 
watching out that it does not intersect the y-parallel lines ^(fc) E^nd Xj^+i), that 
delimit its corresponding subsection. Indeed, such lines can simply be parallel 
to ^(fc) and X(_k+i)- We will introduce them as: 

fik) = ^ , (22) 

where Xk is the common x-coordinate of all points belonging to ^(fe). In the 
same sense, by Xk we will represent the common a;-coordinate of the lines ^P(k)- 
The next step is to collect a finite set of values J over each line V'(fc)- For 
every crossing line, the quantity J of collected values must be large enough to 
warrant that the interpolating functions {/(fc)} (in our example interpolated by 
piecewise-defined polynomials of degree 1) will adequately approach the remain- 
ing conductivity values defined over the line f(k)- This step of the procedure is 



illustrated in Figure 1(b) 



Let us assume that the conductivity inside every subsection is represented 
according to the expression: 

(23) 

Xk + ^(k) 

where Xk denotes the common z-coordinate of the function (p(k)^ introduced in 



(22), f(k){y) is the interpolating function approached with the J values of the 
conductivity, recollected over the Lp(^k) line, and A(^k) is a positive real constant 
such that X + A(k) 7^ within every subsection. 

From this point of view, and supposing that we have K subsections, the 
conductivity a inside the bounded domain f2, can be approached by means of 
the piecewise function: 
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' pw 



( 2^ + ^(1) J- I \ 



a; G [a;2,a;3) 



(24) 



Here, as appointed in Figure 1(a) xq represents the common x-coordinatc of the 
first y-axis parallel line ^(o) employed for subdividing 17, whereas xk represents 
the common coordinate of ^(^k) ■ K is evident that the piecewise function ( 24 1 
is separable- variables. 

Then, according to the Theorem [2j and the relations (17), it follows that 

(lSS-/(2)(y))' : ^e[^i,X2); 



XK+A, 



whereas 



Go 



f{K)iy) 



Xi+-4(i) 


/(I) (a) J 




. 1 ^ 


X(2}+A(2) 


f(2){y) ) 




. 1 ^ 




f{K)(y) J 



X e [xk-i,xk]\ 

X e [a;o,a:i); 
X e [a;i, 2:2); 

X e [xk-i,xk]- 



For the generating pair (i^i, Gi) we will simply have that 

Fi — -^(Tpw, Gi ^ i (-^(Tpw) ; 

These are the generating pairs (Fo, Go) and {Fi, Gi) that we will employ for 
the numerical approach the formal powers. 



3 Numerical solutions of the forward Dirichlet 
boundary value problem. 

In order to evaluate the effectiveness of the piecewise function (Jpw, introduced 



in (24), let us approach the solution of the forward Dirichlet boundary value 
problem corresponding to ([T]), at the perimeter of the unit circle, imposing an 
exact solution as the boundary condition Uc|r- 

Proposition 1 Let 

^Tol) ■ ^^^^ 
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Then the function 



u= +0.1{x + y), (26) 



3 

will be a particular solution of |Ip. 

3.1 Numerical approach of the formal powers. 



We will study a numerical method for approaching elements of the set ( 20 ) . A 
more detailed description of this method, including a variety of special examples, 
can be found in In the next paragraphs, we will focus our attention in the 
construction of the subset of formal powers 



|4")(l,0;z),z(")(l,0;z)r 

^ J n— 



because not any significant alteration is needed when approaching the formal 
powers with coefficient = i- 



Taking into account that the integral operators introduced in (131 are path- 
independent i2jt8j, let us consider a radius R of the unit circle with center at 
zo = 0, as the rectifiable curve 77 described in Q. 

We shall consider P + 1 points equidistantly distributed on R, being the first 
r[0] = and the last r[P] = 1: 

Thus we can construct a set of coordinates according to the formulas: 

x[p] = r[p]cose'g, 
y[p] = r[p]sin6',; 

where 9q is the angle corresponding to R. 



According to (17 1, the data (28) will be used to obtain the sets of values 



Fo[z[p\) = {y[p? + 0.1) ; {x[pY + 0.1) % 

F^{z[p]) = {y[pf + Q.ly'^ (xb]^ + O.l)"' ; 

where z[p\ = x[p\ + iy[p]. Their associated functions Go(z[p]) and Gi{z[p\) will 
be constructed according to (|8]), whereas the adjoin pairs {Fq{z[p\),G'^{z[p])) 
and {Fl{z[p\), G\{z[p\)) will have the form introduced in the Definition[3j 
From ([l2|, it immediately follows that 



Z^^\l,Q-z[p]) - Fo(^b]), 
zf\l,Q;z[p]) - Fi(z[p]); 



(29) 



Hereafter, each formal power with n > will be always approached considering 
P + 1 equidistant points within the interval [0, 1]. 
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Consider the formal powers Zq"'\i,0; z[p]). To approach their values at 
the points z[p\ — x[p\ + iy[p\, we shall employ a variation of the trapezoidal 
integration method: 

4"'(l,0;zH) = 

SFoiz[p]) ■ Re J: (zf 0; z[s + 1]) • {z[s + 1])) dz[s] + 

p— 1 

+5Go(z[p]) • Re E (^("-^^(1, 0; z\s + 1]) • Fq* (z[,s + 1])) + 



+,5Go(z[p])Re £ (zf ^^(1, 0; z[,s]) • Fq* dz[s]; 

where 

d2:[s] = (z[s + 1] - z[s]) , 

and (5 is a real constant factor, heuristically selected, that contributes to the 
numerical stability of the method. In general, every experiment requires an 
individual set of trials in order to estimate an adequate value of b. Particularly, 
for the examples further considered in this work, (5 = 9 provided de best results. 

It is also important to remark that once we have adopted the expression 
(30) for approaching the formal powers, we implicitly employ a piecewise inter- 
polating polynomial function of degree 1, to relate every value Z^"^(l,0;z[p]), 
for 59 = 0, 1, P; and n = 0, 1, N. 

Notice also that, according to the third property of the Theorem [3] for 
Vn > 0, we have that 

Z(")(l,0;z[0])^0. 

Iterating the last procedure we can approach + 1 formal powers that will 
conform the set 

|rcZ("^(1,0;z)|'^ ; 

and after making minor modifications, we shall also approach the set of N formal 
powers 

|ReZ^")(i,0;z)|'^ ; 

pointing out that, according to the Theorem [sj IieZ^^\i,0; Zr) = 0. 

Performing the full procedure for a wide enough quantity Q of radii R, each 
one at some angle 9q: 



2n 



Q-i 

^, = 1-^\ . (31) 

9=0 



we will be able to approach 2N + 1 discrete elements of the set ( 20 1 , introduced 
in the Theorem [5| 



fReZ("'(l,0,z)|r|'^ U f ReZ(")(z, 0, z)|r|'^ . (32) 

I. J 71—0 J n—l 
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Finally, executing a classical Gram-Schmidt orthonormalizing process over 
such elements, and applying a standard cubic splines interpolating method, we 
will obtain a set of 2N + 1 orthonormal piecewise continuous functions, defined 
at every point of the boundary F: 

{un{l):leT}ll,. (33) 



Obtention of the constant coefficients for approaching 
the boundary condition. 



3.2 



Once we have build the set ( |33[ ) , there exist a variety of techniques for approach- 
ing the coefficients {a„}^^0) that will accomplish the task of asymptotically at- 
tain the boundary condition Mc|r (see e.g. [7]). Here, we will allocate as many 
points at the boundary as the number of functions contained in ( 33 ) . 

More precisely, we will trace 2N + 1 radii from the center zq = 0, at the 
angles 



27r 



2^^ + lJ„=o 



2JV 



(34) 

where 



and we will evaluate the boundary condition udr on the points (a;„, y„ 
the radii intersect the boundary F, obtaining a set of values {7n}^i^o- 

We will assume that the number N of formal powers is big enough to warrant 
that, given a number e > 0, there exist a set of real coefficients {a„}^^Q such 
that 



j \ ^ anUn(0 - "c|r 



dl 



<e, leT. 



Thus, the following matrix relation must hold 



where 



U 



and 



uoixQ,yo) 
uo{xi,yi) 
uo{x2,y2) 



ui{xQ,yo) 
ui{xi,yi) 

Ul{x2,y2) 



U2N{xo,yo) 
U2N{xi,yi) 

U2N{x2,y2) 



Uo{x2N,y2N) Ui{x2N,y2N) ■■■ U2N{x2N,y2N) 

7^ = [ao;ai;a2; ■ ■ ■ ;a2N] ; 

^ = [7o;7i;72; • • • ■,12n] ■ 



(35) 



for (35) such that 



Given the orthonormality of the set (33), there will exist an inverse matrix 
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3.3 Comparison between the solutions obtained employ- 
ing the original conductivity, and the piecewise ap- 
proached conductivity. 



We will perform the method described in the Subsection 3.1 considering the 
domain il as the unit circle, and a total error £ of the form: 



£ = 



y a„u„(0 - Uc|r^ dlj ■ (36) 



Beside, we shall fix the following parameters: 

1. Maximum number of formal powers N = 30. 

2. Total number of radii Q — 1000, located at the set of angles defined in 
([31). 

3. Total number of points per radius, where the formal powers are defined, 
P = 1000. 

Therefore, we will obtain a set of 61 orthonormal functions defined at the 
boundary T. 

Valuing the boundary condition (26) at the intersecting points of the Q 
radii R with F, and employing the conductivity a defined in (25), the procedure 
detailed in the Subsection 13. 11 reached a total error: 



£i = 3.6786 X 10"^. (37) 

The same procedure was executed considering the piecewise separable- variables 
conductivity function CTpw(2;,2/), introduced in (24 1, with the same parameters 
indicated in Subsection |2.4[ and assuming also that: 



1. The number of subsections K = 1000; and 

2. The constants 

Ak =60, fc = 1,2, ...,1000. 
The resultant total error was: 

£ = 7.7263 X 10"^ 

The magnitude of this total error indicates that (7p^f,{x,y), introduced in 
(24), can effectively approach separable- variables functions, and in consequence, 
we can use it to apply the elements of the Pseudoanalytic Function Theory for 
analysing the forward Dirichlet boundary value problem of ([T]), in the cases 
when the conductivities a do not originally posses a separable-variables form. 
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Conjecture 1 Let a he a conductivity function, defined within a hounded do- 
main VL (M^) , such that: 

a : n (M^) ^ M; 

and let T he the boundary of such domain. Then it is possible to approach a by 
means of a piecewise separable-variables function of the form (24-): 



«{x,y) = < 



a:+A(2) 



f{i){y) 
f{2){y) 



X e [x^),xi); 

X e [XI,X2)] 



f(K){v) ■ X e [xk~i,xk]; 



which can he employed for constructing a finite set of formal powers of the form 



20), in order to approach solutions of the forward Dirichlet boundary value 
problem corresponding to the equation |Ip, when a boundary condition Uc\t is 
imposed. 

The postulate remains a conjecture because not any formal extension of the 
proof posed in [B], about the completeness of the set (20), is known for the case 
of piecewise separable-variables functions within bounded domains. 



3.4 A special case of piecewise separable-variables con- 
ductivity functions. 

It is natural to inquire about the behaviour of the interpolating method de- 



scribed in the Subsection 2.4 when changing the values of the constant param- 
eters K, ^(fc) or J . In this direction, the work [12] provides some information 
about such specific questions. Nevertheless, the validity of the Conjecture [T] 
allows us to study one property that will significantly simplify our analysis. 

Proposition 2 Let Vl (M^) be a bounded domain, and let us denote by T its 
boundary. Every conductivity function 

can be considered the limiting case of a piecewise separable-variables conductivity 



function CTp^ of the form (24-), at every point (x, y) G £7, when the number K of 



over defined in (22), tends to infinity: 



subdomains introduced in Suhsection\2.4\ and the number J of collected values 



a{x,y)^ lim ap^{x,y). (38) 

-K, J— >-oo 



Moreover, since 



lim 



X 



A 



(k) 



K^oo Xk + ^(fe) 

from the Theorem^ it follows that the generating sequence of this limiting case 
will be periodic, with period c = 1 . 
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Proof 1 Let xq be the minimum of the subset of x-coordinates corresponding 
to the points {x, y) e Vt, and let xk be the maximum. We can divide VL into K 
subdomains {^kl^^i, by employing a set of equidistant y-axis parallel lines (see 



Figure l(a)^ : 

Y (k- 1){XK - xo)\^^^ 
X(k-i) ^ xo-\ — 1 . (39) 

Thus, the subdomains will be defined in the following form: 

{nk\x,yeCk (X(fe_i),X(fe))nr!}^^^, 

where Ck represent the set of points {x,y) within the subsection 

of the plane bounded by the pair of y-axis parallel lines and X(j^y 

Let us consider the set of common x-coordinates belonging to the lines ip(^k) 
introduced in {22): 

{Xk}k=i ■ 

From \3M, we have that 



K 



lim |xfe - a;fc_i| = 0; k = l,...,K; 



where Xk are the common x-coordinate of the lines X(fe), and in consequence we 
will have that 

Vx e fife : X ^ Xk] 



It immediately follows 



X 



A, 



lim '^4^ = 1, (40) 

K^oc Xk + A(fc) 

thus every subdomain flk will be conformed by the points 

{Qk I {x,y) e ip^k)r\n}. 

Furthermore, since the number J of conductivity values, obtained by evalu- 
ating a along the line ^P(k)j o,^so tends to infinite, not any interpolation method 
will be required for approaching ff^ (y) . It will simply coincide with the values of 
a evaluated at {xk,y)- 

Finally, it follows from (4-0) that 

lim ^^J^^ ■ fk{y) = fkiv)- 
K^oo Xk + Ak 

Thus, according to the Theorem the generating sequence for numerically ap- 
proaching the subset of formal powers: 

|z(")(l,0;z),z(")(z,0;^)r (41) 

^ J n—Q 

will be periodic, with period c — 1. 
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The last proposition indicates that the full procedure described in the Sub- 
section 3.1 can be performed considering: 

Fo = \fo; Gq = i (Vct)"^ , 
(Fo,Go) = (^^i,Gi). 

Then, using the same parameters N, Q and P, shown at the beginning of the 
Section |3.3[ and after performing the full numerical procedure, the total error 
was 

£2 = 4.2458 X 10"^ (42) 



which is indeed bigger than the one obtained in ( 37 ) , where all operations were 
performed strictly following the postulates of the Pseudoanalytic Function The- 
ory [2]. Nevertheless, the error (42) is small enough for considering we have 
obtained an acceptable approach of the boundary condition. 

Remark 2 Executing the same logical steps described in Subsection \2.4\ but 
dividing the domain fl by a set of x -axis parallel lines 

{y(o),i(i),...,r(K)}, 

we can approach a piecewise separable-variables conductivity function of the form 



<7pw{x, y) 



= < 



V e [yo.yi); 
y G [2/1,2/2); 



V+A(K) 



f{K){x) ■■ y e [yK-i,yK]; 



whose limiting case, according to the Proposition^ will reach a periodic gener- 
ating sequence with period c — 2. More precisely, the Bers generating pairs will 
have the form 

Fo = ^, Go = ^(^/^)-^ 
Fi = iV^)-\ Gi=iV^. 

Performing the full numerical calculations, and considering the same values N , 
Q and P of the Subsection \3.3[ we obtained a total error 

£ ^ 5.0863 X 10"^ 



which is slightly different that the error ( 57), obtained for the limiting case where 
c = 1 . Nevertheless, it is possible to appreciate that the computational resources 
required for analysing the forward problem when c = 2, are bigger than such 
required for c = 1 . Therefore, hereafter we will exclusively utilize the approach 
where c = 1 . 



17 



4 Examples of conductivity functions with exact 
representations. 

We will analyse a selected set of conductivity functions, for which an exact 
solution is known, in order to impose it as the boundary condition wdr to be 
approached. Once more, the experiments are performed within the unit circle. 
Notice that none of these examples posses a separable- variables form. 

4.1 The exponential case. 

Proposition 3 Let us consider the conductivity function 

a = e"^^, (43) 
where a is a real constant. A particular solution of ^ is 

u = e~"^^ (44) 

The Table 1 contains a condensed relation of the errors £ when changing the 
number of total points per radius P, number of radii Q and maximum number of 
formal powers N, considering a — 1. We shall notice that the behaviour of the 
total error does not keep a clear pattern of change, when the other parameters 
are modified. More precisely, the total error does not decrease monotonically 
when the number of points per radius P, or the number of radii Q, do it. 

On the other hand, the total error does decrease when the number N of 
employed formal powers does. Indeed, the magnitudes of the total error indicate 
that this technique is appropriate for solving this boundary value problem. 

The Table 2 contains the errors £ considering a — 5. Their behaviour is, in 
general, similar to the one of the prior Table. 

4.2 The polynomial case. 

Proposition 4 Let us assume that the conductivity function has the form 

a = a{x + y) + 10, (45) 

thus the function 

w = In {a{x + y) + 10) , (46) 

will be a solution of 

Among the other cases, this is a singular example, because utilizing a small 
number of points per radius P, and a small number of radii Q, together with 
a relatively small number of formal powers N , the convergence of the method 
is acceptable. We expose the results obtained for the case when a = 1 in the 
Table 3. 

When considering the case a — 5, the results presented in Table 4 preserve 
the properties remarked in the previous Table, even when only 5 formal powers 
(11 base functions) are used. 
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Mumber of formal powers 

.V 


Number of radii 

P 


Number of points per radius 

Q 


Tot.-,! error 


30 


1000 


1000 


1.9492 X 10"** 


30 


1000 


800 


2.1979 X 10-« 


30 


1000 


600 


2.2281 X 10"** 


30 


1000 


400 


2.2156 X 10-« 


30 


1000 


200 


1.7241 X 10-« 


30 


800 


1000 


1.5221 X 10"** 


30 


600 


1000 


1.7483 X 10-* 


30 


400 


1000 


1.1651 X 10-« 


30 


200 


1000 


2.6761 X 10"** 


20 


1000 


1000 


3.2741 X 10"* 


10 


1000 


1000 


1.9312 X lO-'^ 


30 


500 


500 


1.1009 X IQ-^ 


20 


500 


500 


1.0376 X 10"** 


10 


500 


500 


1.7330 X 10-'^ 


30 


100 


100 


6.6772 X lO-'^ 


20 


100 


100 


6.9343 X 10"^ 


10 


100 


100 


8.3178 X 10-7 


10 


100 


50 


7.8103 X 10-7 


10 


50 


50 


8.8030 X 10-** 


5 


50 


50 


0.0317 



Table 1: Table of values corresponding to the non separable- variables exponen- 
tial case a = e^^. 
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Mumber of formal powers 

.V 


Number of radii 

P 


Number of points per radius 

Q 


Tot.-,! error 


30 


1000 


1000 


3.3167 X lO"'^ 


30 


1000 


800 


3.4754 X 10"^ 


30 


1000 


600 


3.0912 X 10"^ 


30 


1000 


400 


3.3658 X 10-'^ 


30 


1000 


200 


3.3271 X 10"^ 


30 


800 


1000 


2.6301 X 10"^ 


30 


600 


1000 


2.2022 X lO-'^ 


30 


400 


1000 


5.7358 X 10-^ 


30 


200 


1000 


6.4704 X 10"^ 


20 


1000 


1000 


1.6141 X 10-^ 


10 


1000 


1000 


0.1511 


30 


500 


500 


3.5765 X 10^^ 


20 


500 


500 


1.1817 X 10"^ 


10 


500 


500 


0.1067 


30 


100 


100 


7.2363 X 10^5 


20 


100 


100 


1.1286 X 10"* 


10 


100 


100 


0.0450 


10 


100 


50 


0.0210 


10 


50 


50 


0.0261 


5 


50 


50 


9.4212 



Table 2: Table of values corresponding to the non separable- variables exponen- 
tial case a = e^^^. 
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A'' 


P 


Q 


E 


30 


1000 


1000 


3.6530 X 10"^ 


30 


1000 


800 


3.6528 X 10"** 


30 


1000 


600 


3.6515 X 10"** 


30 


1000 


400 


3.6482 X 10"*^ 


30 


1000 


200 


3.6572 X 10"** 


30 


800 


1000 


4.3271 X 10-8 


30 


600 


1000 


3.5882 X 10"** 


30 


400 


1000 


2.1136 X 10"** 


30 


200 


1000 


1.1306 X 10-** 


20 


1000 


1000 


3.2499 X 10"** 


10 


1000 


1000 


5.9790 X 10-** 


30 


500 


500 


2.0991 X 10"** 


20 


500 


500 


3.1329 X 10-« 


10 


500 


500 


4.4110 X 10"** 


30 


100 


100 


2.8376 X 10"** 


20 


100 


100 


6.5198 X 10-« 


10 


100 


100 


1.1667 X lO"'^ 


10 


100 


50 


1.1973 X 10"^ 


10 


50 


50 


1.1999 X lO"'^ 


5 


50 


50 


1.6714 X lO"'^ 


5 


15 


15 


4.0953 X 10"^ 



Table 3: Table of values corresponding to the non separable-variables polyno- 
mial case (J = X + y + IQ. 
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A'' 


P 


Q 


E 


30 


1000 


1000 


1.5315 X lO"'^ 


30 


1000 


800 


1.5330 X 10-^ 


30 


1000 


600 


1.6351 X 10"^ 


30 


1000 


400 


1.5171 X lO"'^ 


30 


1000 


200 


1.5303 X 10- ^ 


30 


800 


1000 


1.0403 X lO"'^ 


30 


600 


1000 


9.3654 X 10"** 


30 


400 


1000 


5.1980 X 10"** 


30 


200 


1000 


3.3674 X 10-** 


20 


1000 


1000 


2.3294 X 10"^ 


10 


1000 


1000 


2.8096 X 10-4 


30 


500 


500 


6.6396 X 10-** 


20 


500 


500 


9.4135 X 10-** 


10 


500 


500 


1.9824 X lO-'^ 


30 


100 


100 


2.6497 X 10-^ 


20 


100 


100 


4.1895 X 10-^ 


10 


100 


100 


8.3724 X 10-^ 


10 


100 


50 


3.1943 X 10-^ 


10 


50 


50 


4.7861 X 10-5 


5 


50 


50 


0.0145 


5 


15 


15 


0.0039 



Table 4: Table of values corresponding to the non separable-variables polyno- 
mial case 0- = 5(x + y) -\- 10. 
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Number of formal powers 

N 


Number of radii 

P 


Number of points per radius 

Q 


Total error 

£ 


30 


1000 


1000 


1.1213 X 10"** 


30 


1000 


800 


1.9372 X 10-* 


30 


1000 


600 


1.7023 X 10-« 


30 


1000 


400 


1.9995 X 10"* 


30 


1000 


200 


1.9013 X 10"** 


30 


800 


1000 


1.6567 X 10-8 


30 


600 


1000 


1.8385 X 10-« 


30 


400 


1000 


1.2829 X 10-8 


30 


200 


1000 


3.7494 X 10-8 


20 


1000 


1000 


1.1449 X 10-^ 


10 


1000 


1000 


4.7292 X 10-'* 


30 


500 


500 


1.3655 X 10^8 


20 


500 


500 


7.4664 X 10-8 


10 


500 


500 


5.5399 X 10--* 


30 


100 


100 


5.9646 X 10-'^ 


20 


100 


100 


6.8813 X 10-^ 


10 


100 


100 


4.5270 X 10-'' 


10 


100 


50 


4.5305 X 10--* 


10 


50 


50 


4.1998 X 10-"* 


5 


50 


50 


0.0125 



Table 5: Table of values corresponding to the non separable- variables Lorentzial 



case a — [{x + y)^ + 1 



4.3 The Lorentzian case. 

Proposition 5 Let the conductivity function have the form 

1 



[x + y) + a 



An exact solution for the equation |7p 

{x + y) 



IS 
3 



3 +a{x + y). (48) 

We will analyse the case when the constant a = 1, showing the obtained 
results in the Table 5. Once more, we detect that decrementing the number of 
points per radius P, as well the number of radii Q, does not necessarily provoke 
a diminution of the total error £. 



The case when a = 0.01, in (47) and (48), does provide additional informa- 
tion for our study. The maximum value of the conductivity cr, at the origin, is 
ten thousands times bigger than the one belonging to the cases studied before. 
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Number of formal powers 

N 


Number of radii 

P 


Number of points per radius 

Q 


Total error 

£ 


30 


1000 


1000 


0.1671 


30 


1000 


800 


0.1671 


30 


1000 


600 


0.1671 


30 


1000 


400 


0.1671 


30 


1000 


200 


0.1672 


30 


800 


1000 


0.1491 


30 


600 


1000 


0.1283 


30 


400 


1000 


0.1030 


30 


200 


1000 


0.0662 


20 


1000 


1000 


0.6420 


10 


1000 


1000 


3.1839 


30 


500 


500 


0.1165 


20 


500 


500 


0.4509 


10 


500 


500 


2.2466 


30 


100 


100 


0.0265 


20 


100 


100 


0.1559 


10 


100 


100 


0.9435 


10 


100 


50 


0.0870 


10 


50 


50 


0.5879 


5 


50 


50 


120.3691 



Table 6: Table of values corresponding to the non separable-variables 



Lorentzian case cr = [{^ + 2/)^ + 0.01 



From this point of view, the notorious increment of the total errors is justi- 
fied. Yet, it is interesting that the diminution of P and Q do not influence that 
much in the errors. The Table 6 contains the numerical data that sustain our 
affirmations. 

4.4 The sinusoidal case. 

This example has been selected because it proved to be one of the most chal- 
lenging cases for testing the performance of the numerical method. Our analysis 
will be based in the following proposition. 

Proposition 6 Let us consider the sinusoidal conductivity 

(7 = 1 + sinaxy. (49) 

We can verify by direct substitution that the function 

«=(tan(^)+l)"\ (50) 

is a particular solution of tw. 
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Number of formal powers 

N 


Number of radii 

P 


Number of points per radius 

Q 


Total error 

£ 


30 


1000 


1000 


1.2451 X 10-« 


30 


1000 


800 


1.2263 X 10-^ 


30 


1000 


600 


1.1838 X IQ-^ 


30 


1000 


400 


1.3812 X 10^8 


30 


1000 


200 


1.4932 X 10"** 


30 


800 


1000 


9.4572 X 10-9 


30 


600 


1000 


1.1218 X IQ-^ 


30 


400 


1000 


9.7996 X 10-9 


30 


200 


1000 


9.6253 X 10-9 


20 


1000 


1000 


2.3682 X 10-^ 


10 


1000 


1000 


3.9153 X 10-^ 


30 


500 


500 


9.7774 X 10^9 


20 


500 


500 


2.4493 X IQ-^ 


10 


500 


500 


5.6988 X 10-^ 


30 


100 


100 


3.7062 X 10-^ 


20 


100 


100 


5.1172 X 10-^ 


10 


100 


100 


9.3918 X 10-'^ 


10 


100 


50 


9.3414 X 10"^ 


10 


50 


50 


6.9851 X 10-5 


5 


50 


50 


0.0497 



Table 7: Table of values corresponding to the non separable- variables sinusoidal 



case a — 1 + sinxy. 



Let us focus our attention into the case a = 1. The Table 7 contains what 
it could be considered a set of acceptable approaches for the boundary value 
problem, displaying a similar behaviour to the rest of cases previously shown. 
More interesting information reaches when considering a = 5 in the argument 



of the sinusoidal expression ( 49 ) , but keeping the value a = 1 in the boundary 
condition (50). This is mainly because the condition is not defined at several 
points of the boundary F, when considering a = 5. 

Even not any formal extension of the Theorem [5] is known for the limiting 
case posed in the Proposition [2j about the completeness of the set (20), one 
could expect a diminution of the total error £ when increasing the number N of 
formal powers, as registered in the other examples. But the Table 8 shows that, 
for this particular case, the behaviour is the opposite. This implies that we have 
found an example that could be very significant to adequately understand the 
properties of the formal powers constructed by virtue of the Proposition [2] The 
case must certainly be analysed with more detail in further works. 

Finally, we shall remark that 

a — 1 + sin 5x2/ > 0: '^^^ U ^ ^- (51) 
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Number of formal powers 

N 


Number of radii 

P 


Number of points per radius 

Q 


Total error 

£ 


30 


1000 


1000 


1.0338 X lO'^ 


30 


1000 


800 


1.2828 X 10^ 


30 


1000 


600 


7.1245 X 10^ 


30 


1000 


400 


2.7315 X lO'^ 


30 


1000 


200 


1.4515 X 10* 


30 


800 


1000 


2.7273 X 10^ 


30 


600 


1000 


7.1022 X 10^ 


30 


400 


1000 


2.2060 X 10^ 


30 


200 


1000 


1.9391 X 10^ 


20 


1000 


1000 


3.6205 X 10"^ 


10 


1000 


1000 


2.7627 X 10^ 


30 


500 


500 


2.3643 X 10^ 


20 


500 


500 


1.2936 X 10"^ 


10 


500 


500 


7.0407 X lO'^ 


30 


100 


100 


3.2401 


20 


100 


100 


116.1873 


10 


100 


100 


2.2541 X 10^ 


10 


100 


50 


0.1006 


10 


50 


50 


24.2426 


5 


50 


50 


47.5417 



Table 8: Table of values corresponding to the non separable- variables sinusoidal 



case a = 1 + sin bxy. 



5 Conductivities corresponding to geometrical 
distributions. 

This section is dedicated to analyse what it could well be considered one of the 
most important contributions of this work to the State of the Art for solving 
the Dirichlct boundary value problem of ([T]), employing the elements of the 
Pseudoanalytic Function Theory: The study of conductivity distributions rising 
from geometrical cases. 

As it was mentioned before, the tools provided in ^ and ^8] impose an 
important condition, seldom fulfilled in physical experimental models: The con- 
ductivity a must be a separable- variables function. The Proposition [2] allows us 
to overpass this restriction by stating that any conductivity function: 

a : n (M^) ^ M, (52) 

can be considered the limiting case of a piecewise-defined function, of the form 
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( 24 ) . The Section above studied the cases when the a possesses an exact repre- 
sentation, but the mapping ( 52 ) indicates that it is enough to know the value of 
a at every point {x^y) within the domain f2. Thus, the conductivities given by 
geometrical distributions, can be treated with the same tools that we employed 
in the Section |4l 

Before continuing with this branch of the analysis, we shall remark that, as 
it was explained in Subsection 3.1 when employing the expression (30) for nu- 
merically approaching the formal powers, we are indeed performing a piecewise 
interpolating polynomial method of degree 1, on every pair of valued-points on 
the radius where the conductivity a is defined. 

This implies that the discontinuities of the first kind, appearing in the con- 
ductivity over one radius, will not be taken into account, since there will always 
exist a line with finite slope, relating the values a[r[p\) and a{r[p + 1]), between 
which the discontinuity is found. 

Nevertheless, the set of examples presented in the upcoming Subsections, 
will show that this approach provides numerical solutions with an acceptable 
degree of accuracy. 



5.1 The concentric disks. 



Let us consider a piecewise conductivity function, in polar coordinates, of the 
form: 



(53) 



Here r denotes the radius. The Figure 2 illustrates this case. To select a 
boundary condition for the conductivities defined by geometrical figures, with- 
out performing physical measurements, it is not a trivial task. Nevertheless, 



100 


re [0,0.2); 


30 


r e [0.2,0.4) 


20 


r e [0.4,0.6) 


15 


r e [0.6,0.8) 


10 


r e [0.8,1]. 



the conductivity defined in ( 53 ) may be related with the separable- variables 
Lorentzian case, studied in Proposition [l] For this reason, the boundary condi- 
tion will be a variation of the expression (|26|) : 



O.Ol{x + y). 



The Table 9 summarizes the total errors £ obtained when changing the 
number of points P, the total radii Q, and the employed formal powers N. The 
behaviour of the total error is very similar to the one exhibited in the Lorentzian 
case of the Subsection 14.31 



5.2 Disk out of the center, within the unit circle. 

Let us consider now the conductivity function illustrated in the Figure 3. The 
blue section represents conductivity values of a = 10, whereas the red area 
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Figure 2: Conductivity defined according to the expression (53) 



Number of formal powers 

N 


Number of radii 

P 


Number of points per radius 

Q 


Total error 

£ 


40 


1000 


1000 


3.6234 X 10-'' 


40 


1000 


800 


2.7852 X 10"'-' 


40 


1000 


600 


2.5213 X IQ-'-' 


40 


1000 


400 


2.3199 X IQ-^ 


40 


1000 


200 


1.4331 X 10-^ 


40 


800 


1000 


3.6234 X IQ-'-' 


40 


600 


1000 


3.6234 X IQ-'-' 


40 


400 


1000 


3.6234 X 10-y 


40 


200 


1000 


3.6234 X 10-9 


20 


1000 


1000 


3.3615 X 10"'-* 


20 


500 


500 


2.2764 X 10-9 


40 


100 


100 


2.4841 X 10-^ 


20 


100 


100 


3.3996 X 10-7 


10 


50 


50 


8.2633 X 10-*^ 


5 


50 


50 


1.1721 X 10-^ 



Table 9: Table of values corresponding to the conductivity a of the Figure 2: 
The concentric disks. 
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Figure 3: Conductivity function with one disk within the unit circle. The red 
surface represent a — 100, whereas the blue sections denote a = 10. 




-1 -0.5 0.5 1 



Figure 4: Square within the unit circle. The blue surface represents a conduc- 
tivity (7 = 10, and the red indicates a = 100. 



expresses a =^ 100. The red disk corresponds to the equation 

(x- 0.6)2 +2/2 < Q 2. 

For this case, and considering once more the possible similarity with the 



Lorentzian case ( 26 ) , we will impose the boundary condition: 



{x - o.ey + 



+ 0.01 (x - 0.6 + y) . 



The results of the numerical calculations are displayed in the Table 10. It 
is remarkable that the total errors are, in general, five thousands times bigger 
than such shown in the Table 9, where the concentric disks were analysed. 



5.3 Square within the unit circle. 

As plotted in Figure 4, inside the unitary circle we located a perfect square, 
whose superior and inferior sides are parallel to the x-axis, whereas its left and 
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Number of formal powers 

N 


Number of radii 

P 


Number of points per radius 

Q 


Total error 

£ 


40 


1000 


1000 


7.8082 X lO-'^ 


40 


1000 


800 


6.9845 X 10-* 


40 


1000 


600 


5.9861 X lO"'^ 


40 


1000 


400 


4.3758 X 10-'' 


40 


1000 


200 


3.1089 X 10-* 


40 


800 


1000 


7.8048 X 10-* 


40 


600 


1000 


7.7999 X 10-* 


40 


400 


1000 


7.6313 X 10-* 


40 


200 


1000 


7.8686 X 10-* 


20 


1000 


1000 


1.6829 X 10-3 


20 


500 


500 


1.8330 X 10-3 


40 


100 


100 


6.1065 X 10-^ 


20 


100 


100 


7.7260 X 10-* 


10 


50 


50 


1.5540 X 10-3 


5 


50 


50 


3.9329 X 10-3 



Table 10: Table of values corresponding to the conductivity a of the Figure 3: 
One disk out of the center of the unit circle. 



right sides are parallel to the y-axis, being their lengths all equal to 0.65. The 
square corners are equidistant to the origin, and we have fixed one point of the 



set (27) at each corner, so they can effectively take part into the calculations of 
the formal powers. For this case, we will impose a boundary condition of the 
form: 

3 3 

u=^^+0.1(x + y). (54) 

The obtained errors are shown in Table 11, and even they are, in general, 
one thousand times bigger than the errors calculated for the disk out of the 
center, presented in the previous Subsection, it is important to remark that the 
presence of the corners did not influence significantly the convergence of the 
approach on the boundary T. 

This implies that, for this conductivity distribution, with the boundary con- 



dition ( 54 ) , not any supplementary regularization method is required. 

As a matter of fact, the complete following Section is composed by a col- 
lection of cases where not any additional regularization methods are needed to 
successfully approach the boundary condition, in non-smooth bounded domains. 



6 Brief analysis of a non-circular domains. 

Many numerical methods, for solving partial differential equations at some do- 
main 57 in the plane, often require regularization techniques when the derivative 
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Number of formal powers 

N 


Number of radii 

P 


Number of points per radius 

Q 


Total error 

£ 


40 


1000 


1000 


1.4598 X 10"^ 


40 


1000 


800 


1.4603 X 10^2 


40 


1000 


600 


1.4601 X 10-^ 


40 


1000 


400 


1.4619 X 10-2 


40 


1000 


200 


1.4638 X 10-2 


40 


800 


1000 


1.4160 X 10-2 


40 


600 


1000 


4.0513 X 10^2 


40 


400 


1000 


2.7287 X 10-2 


40 


200 


1000 


1.4638 X 10^2 


20 


1000 


1000 


2.3516 X 10^2 


20 


500 


500 


1.6912 X 10-2 


40 


100 


100 


3.2869 X 10-^ 


20 


100 


100 


3.4552 X 10-2 


10 


50 


50 


6.8926 X 10-2 


5 


50 


50 


1.9362 X 10-1 



Table 11: Table of values corresponding to the conductivity a of the Figure 4: 
A square within the unit circle. 



of the parametric curve describing F, possesses discontinuities (see e.g. [TO]). 
In general, the method posed in this work is not the exception. Nevertheless, 
there exist some particular cases such that not any regularization procedures 
are required, for adequately approaching the boundary condition in non-smooth 
domains. 

The examples we present in the forward paragraphs were heuristically found, 
and it is important to clarify that the convergence of the approached solutions 
was only achieved when establishing a certain set of values. More precisely, we 
fixed the number of radii Q — 100, whereas the number of points per radius 
P = 100. 

It is interesting to remark that, if increasing the number Q in the selected 
examples, the convergence at the non-smooth regions will not be achieved. Nei- 
ther will do if decrementing it. Moreover, we found out that the number of 
formal powers 2N -\- 1 employed in the approach must be equal to 91 (46 formal 
powers with coefficient a„ = 1 and 45 with coefficient a„ = i), for obtaining the 
minimum total error £. This will be valid for all examples shown below. 

This Section is included mainly for remarking that, the numerical solving of 
Dirichlet boundary value problems for elliptic equations on the plane, employing 
the elements of the Pseudoanalytic Function Theory, possesses the property of 
approaching adequately the boundary condition in non-smooth domains, when 
certain conditions are fulfilled. This motivate a deeper study of such property, 
for correctly understanding which are the conditions to warrant the convergence, 
and for which cases this can be effectively done. 



31 



-1 -0.5 0.5 1 1.5 
x-axis 

Figure 5: Example of a conductivity a — (a;^ + O.l) ^ (jj^ + O.l) ^, whithin a 
non-circular domain. 



We shall analyse the behaviour of the method on the boundary F of a non- 
circular domain, that can be described as follows: 

1. The domain coincides with the unit circle at the left-hand side of the 
y-axis parallel line x — cos . 

2. At the right-hand side of this axis, the domain is the area bounded by 
the line segments: 

y = -0.5629a; + 0.8443, , . 

y = 0.5629a; + 0.8443; ^ ' 

in the closed interval x G [cos j^, 1.5] . 
6.1 The Lorentzian case. 

Let us consider a conductivity a within the domain described in the Subsection 
[6] with the form: 

(56) 



0.1 J + 0.1 

imposing as the boundary condition the exact solution: 

u^ ^ > +0.l{x + y). (57) 

The illustration of this example is displayed in Figure 5. 

As mentioned above, we shall fix the number of radii Q — 100, and the 
number of points per radius P = 100. The Figure 6 displays the approach of 
the boundary condition (57 1 using 2A^+ 1 = 91 base functions. We remark that 
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1.5 




-2-1012 
Perimeter of the non-circular domain 
broaclned in anguiar perspective witlnin tine segment 
[-3.1416,3.1416] 



Figure 6: Comparison between the boundary condition (571 and the numerical 
approach employing 2N + 1 = 91 base functions, for the case when a possesses 
the form a = {x^ + O.l)"^ {y^ + 0.l)~\ 



the Figure 6 was traced keeping an angular perspective, broaching an angle 
9 e [— TT, tt], represented by the ordinate-axis, whereas the abscissa-axis plots 
the electric potential values related to each angle. The total error £ for this 
example was: 

£ = 1.8355 X lO"''. 



We also mention that one of the points where the boundary condition (571 
has been valued, was located precisely at the intersection of the line segments 
defined in ( |55[ ): (1.5,0), where we find the most representative non-smoothness 
of the boundary F. This point coincides with the radius corresponding to the 
angle w = 0. 

Another two boundary value points were fixed at (1.5 — cos j^, sin ^q) and 
(1.5 — cos — sin yI^), where the derivative of the parametric curve describing 
the boundary F, is also discontinuous. These points are associated with the 
radii located at the angles to = ^ and a; = — . 

When reducing the number of base functions at 2N + 1 = 51, we observe a 
considerable divergence of the approached solution, around the boundary value 
related to the angle a; = 0, as it is displayed in Figure 7. Once more, we remark 
that one boundary value point was fixed at (1.5, 0). On the other hand, since the 
distribution of the boundary value points is performed according to the angular 
expression (34), the points (1.5 — cos f^, sin f^) and (1.5 — cos fg, — sin f^) were 



not included in the analysis of this example. 

The total error employing 2N -1-1 = 51 base functions was: 

£ = 0.3869. 
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1.5 




■^-3 -2-10123 



Perimeter of the non-circular domain 
broaclned in anguiar perspective witlnin tine segment 
[-3.1416,3.1416] 



Figure 7: Comparison between the boundary condition (571 and the numerical 
approach employing 2N + 1 = 51 base functions, for the case when a possesses 
the form a = (x'^ + QAV^ (y^ + 0.l)"\ 



6.2 Concentric circles. 

We shall consider a geometrical variation of the conductivity distribution in- 



troduced in ( 53 1 , within the same domain 11 described at the beginning of the 
Section [6] as shown in Figure 8. For this example, we will perform the initial 
experiment taking into account 2N -|- 1 = 91 base functions. We also fix three 
boundary value points at the discontinuities of the derivative of F indicated 
before. 



We choose to impose the boundary condition employed in (57 1, and the 
Figure 9 illustrates the behaviour of the approached solution when compared 
with the boundary condition. The total error was: 

£ = 3.4217 X 10"''. 

The second experiment took into account only 2N -I- 1 = 51 base functions. 
Again, the boundary value points were distributed at the perimeter of the do- 



main fl according to the angular model (34), empathizing that one point was 
fixed at (1.5,0). As displayed in figure 10, the highest divergence of the ap- 
proached solution is located around this point. 



6.3 A squared figure inside the non-circular domain. 

As a final trial, let us review the case when the conductivity a corresponds to 
the distribution plotted in Figure 11. This could well be considered a significant 
example for studying the effectiveness of the posed numerical method. 

We find three points where the derivative of the parametric curve, describing 
F, is discontinuous, and we have forced one boundary value point to be located 
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-0.5 
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x-axis 
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Figure 8: Example of a conductivity a composed by concentric circles in a 
non-circular domain. 
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Perimeter of the non-circular domain 
broaclned in anguiar perspective witlnin tine interval 
1-3.1416.3.14161 



Figure 9: Comparison between the boundary condition (571 and the numerical 
approach employing 2N -I- 1 = 91 base functions, for the case when a possesses 
the form plotted in Figure 8. 
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-2-1012 

Perimeter of the non-circular domain 
broaclned in anguiar perspective witlnin tine interval 
1-3.1416.3.14161 



Figure 10: Comparison between the boundary condition (571 and the numerical 
approach employing 27V -f 1 = 51 base functions, for the case when a possesses 
the form plotted in Figure 8. 
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-1 -0.5 .0.5 1 1.5 
x-axis 

Figure 11: Example of a conductivity a with a squared figure inside fi. 
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Perimeter of the non-circular domain, 
broaclned in anguiar perspective witlnin tine interval 
[-3.1416,3.1416] 



Figure 12: Comparison between the boundary condition (571 and the numerical 
approach employing 2N + 1 = 71 base functions, for the case when a possesses 
the form plotted in Figure 11. 



at each one of them. Beside, the figure within the domain fl possesses four 
corners. More precisely, it is a perfect square whose characteristics are the 
same that those employed in the Subsection |5.3| The first example considered 
2N -I- 1 = 91 base functions, being its total error: 



£ = 0.0010. 

The comparison between the boundary condition and the approached solu- 
tion is displayed in Figure 12. The last trial was performed utilizing 2N + 1 = 71 
base functions, forcing a boundary value point at (1.5,0). 

In the proximity of this point the divergence of the solution is considerable, 
as displayed in Figure 13, an this is also shown by the increased total error, that 
for this case was: 



£ ^ 2.4867. 



Appointing 1 The numerical methods used along this work were fully devel- 
oped m GNU C/C++ Compiler, employing a CPU INTEL® processor XEON® 
E5620x&ABm4 CHz, on SLACKWARE 13.37 LINUX operating system. 

The experimental procedures showed that the numerical results can variate 
when using different platforms based in 32B and 645 processor architecture, or 
compilers between Microsoft® Windows® and LINUX operating systems. 

If the reader wishes to perform his own numerical trials, please contact the 
authors to obtain the resource codes. 
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